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ABSTRACT 

The SCARE (Structural Ceramics Analysis and Reli- 
ability ^valuation) computer program on statistical 
fast fracture reliability analysis with quadratic ele- 
ments for volume distributed imperfections is enhanced 
to include the use of linear finite elements and the 
capability of designing against concurrent surface flaw 
induced ceramic component failure. The SCARE code is 
presently coupled as a postprocessor to the MSC/NASTRAN 
general purpose, finite element analysis program. The 
improved version now includes the Wei bull and Batdorf 
statistical failure theories for both surface and vol- 
ume flaw based reliability analysis. The program uses 
the two-parameter Weibull fracture strength cumulative 
failure probability distribution model with the prin- 
ciple of independent action for polyaxial stress 
states, and Batdorf's shear-sensitive as well as shear- 
insensitive statistical theories. The shear-sensitive 
surface crack configurations include the Griffith crack 
and Griffith notch geometries, using the total critical 
coplanar strain energy release rate criterion to pre- 
dict mixed-mode fracture. Weibull material parameters 
based on both surface and volume flaw induced fracture 
can also be calculated from modulus of rupture bar 
tests, using the least squares method with known speci- 
men geometry and grouped fracture data. The surface 
flaw reliability prediction uses MSC/NASTRAN stress, 
temperature, and external boundary area output, 
obtained from the use of linear or quadratic shell and 
three-dimensional isoparametric finite elements. The 
statistical fast fracture theories for surface flaw 
induced failure, along with selected input and output 
formats and options, are summarized. An example prob- 
lem to demonstrate various features of the program is 
i ncl uded. 

INTRODUCTION 

Brittle structures characteristically exhibit a 
large variation in fracture stress which must be taken 
into account in design. This variation in observed 
strength of nominally identical components is due to 
the presence of invisible material imperfections. 


These imperfections are assumed to have a distribution 
in strength and the structure is assumed to fail when 
the strength of the weakest flaw or link is exceeded. 
Structural ceramics are known to contain at least two 
concurrent types of material imperfections or flaw 
populations. One type of unavoidable flaws arise from 
material processing, which are usually restricted to 
the interior of the structure and are referred to as 
volume or intrinsic flaws. Another flaw population 
with distinctly different structural response exists 
on the surface of a ceramic component. These external 
surface or extrinsic flaws can arise from qrinding or 
other finishing operations, or maybe the result of 
environmental factors or of the intrinsic porosity 
intersecting the external surface. Statistical fast 
fracture models, based on weakest link theories (WLT) 
of Weibull, have been previously developed for both 
types of flaw populations (Weibull, 1939; Rufin et al., 
1984; Wertz and Heitman, 1980; Batdorf, 1978; Johnson, 
1983). Recently, a public domain general purpose 
reliability code called SCARE has been generated 
(Gyekenyesi, 1986) to predict the fast fracture 
response of ceramic structures due to volume distrib- 
uted flaws. It is the purpose of this paper to 
describe enhancements to the SCARE program which will 
allow reliability analysis with linear or quadratic 
elements of ceramic components due to the presence of 
a concurrent surface distributed flaw population. 
Consequently, fracture caused by either internal or 
external flaws can then be accounted for in the design 
process. 

The first probabilistic method used to account for 
the scatter in fracture strength of brittle materials 
was introduced by Weibull (1939). His model was based 
on the weakest link theory and required certain statis- 
tical parameters to phenomenologically describe the 
failure response of a material. These statistical 
parameters were usually determined from uniaxial ly 
loaded, simple qeometry specimens. To predict mate- 
rial behavior in other stress states using statistical 
parameters from uniaxial tests, Weibull proposed calcu- 
lating the risk of rupture by averaqing the tensile 
normal stress in all directions. Since this approach 
is arbitrary and requires numerical modeling, other 
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approaches were subsequently introduced. The most 
widely used of these is the principle of independent 
action (PIA) model, which is based on the assumption 
that the principal stresses act independently. The PIA 
fracture theory is the statistical version of the maxi- 
mum stress failure theory, which was widely used in 
early metallic structure design. Most ceramic compo- 
nents in the government sponsored advanced gas turbine 
( AGT) programs used either or both of the above 
described polyaxial statistical failure models. Both 
of these fracture theories, however, can lead to unsafe 
estimates of failure (Gyekenyesi, 1986) since they both 
neglect the shear force, and in case of the PIA hypoth- 
esis, the effects of combined principal stresses. 

In order to introduce a mechanistic fracture cri- 
terion into a statistically based analysis model to 
more accurately predict the onset of catastrophic crack 
propagation, Batdorf developed a new theory which also 
applies to brittle materials in which the cracks are 
confined to the surface (Batdorf, 1973; Batdorf and 
Heinisch, 1978). He assumed that the cracks are ran- 
domly oriented, that they are Derpendicul ar to the 
boundary and that they do not interact. The sizes of 
the cracks are not explicitly treated and fracture 
occurs when the remote macroscopic stress normal to the 
crack plane exceeds some critical stress, a cr , charac- 
terizing that particular crack. Since the cracks are 
assumed normal to the surface, crack orientation is 
given by a single angular measure rather than two, as 
in the case of volume distributed flaws. This simpli- 
fication is not always invoked (Shetty et al.,. 1984), 
although there is strong experimental evidence that it 
is justified (Batdorf, 1973). 

In adding the influence of shear loading on the 
crack face, the external flaws are currently modeled 
as Griffith cracks and Griffith notches. Unlike with 
internal flaw analysis in SCARE, only the total criti- 
cal coplanar strain energy release rate, G c , criterion 
is used to predict mixed-mode fracture (Samos, 1982), 
and the maximum tensile stress criterion, as described 
by Batdorf and Heinisch (1978), is not implemented for 
surface cracks in the program. It has also been con- 
cluded (Batdorf and Heinisch, 1978) that the Griffith 
crack is not an appropriate model of a typical surface 
imperfection, and is included here mainly for histori- 
cal and for academic reasons. 

Surface flaw based reliability analysis has been 
implemented in both the SCARE1 and SCARE2 versions of 
the postprocessor program (Gyekenyesi, 1986). SCARE1 
uses only elemental centroidal principal stresses to 
calculate reliability. In the SCARE2 version of the 
code, all linear or quadratic QUAD8 shell elements are 
further discretized into 9 subelements, which are then 
used with interpolated principal stresses to perform 
all analysis. In order to identify external bound- 
aries, calculate their surface areas and obtain corre- 
sponding surface stress states, appropriate MSC/NASTRAN 
shell elements of negligible stiffness are used, 
together with the previously selected three-dimensional 
HEXA and PENTA elements. The selective use of shell 
elements permits the identification of only those 
external areas which contain potentially failure- 
causing flaws, and it can ignore boundaries from which 
fracture is not likely, such as finite element model 
symmetry planes or compressi vely loaded external 
surfaces. 

PROGRAM CAPABILITY AND DESCRIPTION 

The architecture and basic computational elements 
of the postprocessor program are described by 
Gyekenyesi (1986). Figure 1 shows the flowchart for 


the new surface flaw SCARE2 reliability analysis code, 
which is designed to be identical in its sequence of 
calculations to the previously developed volume flaw 
analysis. For computational efficiency, all experi- 
mental fracture stresses as well as all elemental 
principal stresses are normalized. When using the 
Weibull scale parameter, o 0 , as the normalizing 
factor, it is important to note that different values 
of o 0 are generally used in the volume and surface 
flaw calculations. For efficient transfer of 
MSC/NASTRAN outDUt data to SCARE, FORTRAN logical units 
3, 4, and 7 are used. Input to SCARE is handled 
through logical units 3, 5, and 7 while the output is 
stored on unit 1. 

Presently, the proqram permits use of several 
fracture criteria, crack configurations, and tempera- 
ture dependent statistical material parameters. Uni- 
axial, surface flaw induced, fracture data alonq with 
specimen geometry from four-point bend tests can be 
used to calculate Weibull parameters and the Batdorf 
flaw density coefficient. Fiqure 2 contains the avail- 
able options in fracture criteria and flaw configura- 
tions used to model surface imperfections. Note that 
two of the failure criteria are for shear-insensitive 
cracks, even in multidimensional stress states. The 
other two are used for the shear-sensitive models, 
where the mechanistic G c fracture criterion is 
employed to predict component reliability. Amonq the 
available criteria and crack configurations shown in 
Fiq. 2, the Griffith notch with the G c criterion 
has the hiqhest shear-sensitivity and hence, the high- 
est failure probability for a qiven structure, while 
the PIA approach yields consistently the lowest failure 
estimate. It should also be noted that the Batdorf 
shear-insensitive fracture model is conceptually iden- 
tical to the oriqinally proposed Weibull normal stress 
averaging method, althouqh it has been recognized early 
that, when using this method, an amplification factor 
is required to accentuate the materials sensitivity to 
surface defects (Paluszny and Wu, 1977). No such fac- 
tors, however, are used in the current version of the 
SCARE code. 

In selecting thin shell elements to determine 
surface stress states, the existence of only two sur- 
face principal stresses is intrinsic to the analysis. 

In order to use WLT, any compressive elemental princi- 
pal stress is not permitted to exceed three times the 
tensile principal stress in absolute value. Otherwise, 
compressive stress state predominates and the corre- 
sponding reliability is set eaual to unity. Just like 
in uniaxial compressive loadinq, when using the PIA 
model with Weibull statistics, fracture due to com- 
pression is inadmissible. 

INPUT INFORMATION 

The SCARE computer program requires output from 
MSC/NASTRAN elastostatic analysis to determine compo- 
nent fast fracture reliability. Gyekenyesi (1986) 
discusses some of the available rigid format solution 
sequence options and also how the NASTRAN and SCARE 
programs interface on the Lewis laboratory computer 
system. The surface flaw analysis capability uses only 
HEXA, PENTA, QUAD8, and TRIA6 MSC/NASTRAN library ele- 
ments, and model inq with axisymmetric elements is not 
permitted for external flaw reliability predictions. 
This restriction is due to NASTRAN finite element prop- 
erties which rule out mixinq of axisymmetric elements 
with any other element type in a qiven mesh. Only the 
HEXA and PENTA isoparametric elements should contribute 
to structural stiffness and the planar shell elements 
should have negligible thickness, t, specified in 
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NASTRAN (t = 10~ 6 in. for the rotating disk example 
in this paper). In addition to the previously listed 
reasons for using shell elements, modeling a ceramic 
structure with them will also avoid the need for ele- 
mental stress transformations when calculating required 
surface principal stresses. Fur avail able stress 
recovery options, users should consult the appropriate 
analysis code manuals. With today's automatic mesh 
generation programs, little difficulty will be encoun- 
tered in adding the required external elements. A 
solid element face is identified as an external sur- 
face when its nodes are shared by a shell element. The 
midside nodes may or may not be present, but current 
capability requires that either they are all there con- 
sistently, or they are all absent. Mixing of linear or 
quadratic elements in a given finite element mesh is 
not permitted. When using shell elements, only the 
membrane properties should be invoked and the uncoupled 
bending stiffness should not be included in the model. 

SCARE input requirements for external flaw analy- 
sis are almost identical to the input data for internal 
flaw analysis, which has been previously summarized 
(Gyekenyesi, 1986). The postprocessor program is cap- 
able of doing the volume and surface flaw analysis in 
one execution. However, only the surface distributed 
crack capability will be described in this paper. As 
before, program input can be grouped into three cate- 
gories. The first category, shown in Fig. 3, is called 
the master control deck, which includes organization of 
control indices for fracture criteria, finite element 
model requirements, and material parameter format. 

Note that for calculating statistical parameters, data 
only from four-point bend test can currently be used. 

In addition, two-lines or cards of data entry are now 
required rather than one as in the previous volume flaw 
analysis. The fracture criteria and flaw shapes for 
representing volume and surface imperfections can be 
independently specified, and users can select to per- 
form analysis for volume flaw based fracture only, for 
surface flaw based fracture only, or for both simul- 
taneously in one execution. Details on included 
options and limitations can be found in the program 
user's manual . 

The second category, called the specimen deck, 
uses fracture specimen data required in calculating 
statistical fracture parameters, or direct input of 
material properties including Poisson's ratio. These 
material parameters are generally temperature depend- 
ent. Therefore, provision is made for input of multi- 
ple flaw populations at different temperatures. Unless 
material parameters are directly read, this section 
requires four-point bend bar geometry and surface flaw 
based failure stresses of the sample population. These 
extreme fiber failure stresses must be arranged in 
ascending order and currently up to 200 specimens can 
be used for a given temperature. The specimen deck 
also requires test temperatures, arranged in ascending 
order for multiple temperatures, since calculated mate- 
rial parameters are interpolated within SCARE. At a 
specified temperature, fracture stresses must be unique 
and multiple values of identical magnitudes are not 
permitted. The number of available fracture readings 
for all temperature tests must be the same. Additional 
explanation of the required input, including size limi- 
tations, can be found in the user's manual. 

The last SCARE input category, called the struc- 
tures deck, contains results of the finite element 
structural analysis required for failure probability 
predictions. These include solid element volumes, 
shell element areas, nodal temperatures (MSC/NASTRAN 
does not permit access to element temperatures), ele- 


mental principal stresses and corresponding element 
identification numbers. In the present version of 
SCARE, which relies on MSC/NASTRAN output files, all 
of this data is internally manipulated throuqh sub- 
routine ELEM, and the structures deck requires no 
5yGc t • ic input fcy the user. 

OUTPUT INFORMATION 

Grid point temperatures in a finite element mesh 
of a selected component can be obtained at transient 
or steady state conditions from performing MSC/NASTRAN 
thermal analysis. Combining the resultinq thermal 
loads with additional mechanical loads, the displace- 
ment and elastic stress fields can then be calculated. 
For surface flaw based fracture analysis, shell element 
stresses, areas and temperatures are the required out- 
put from the analysis code. For available output 
options, users should consult the appropriate program 
manuals. The shell element stresses as well as the 
nodal temperatures from MSC/NASTRAN are stored in 
specified punch files, while surface element areas are 
listed in the proqram output file. Consequently, 
NASTRAN 's punch option must be invoked to store the 
required stress and temperature data with the program's 
parameter call feature needed to calculate the required 
elemental areas. 

The first part of all SCARE output data contains 
an echo of important NASTRAN finite element analysis 
results. Identifying labels, element type, and number 
of elements in the model are noted. Since in a larqe 
finite element mesh the stress output could be exces- 
sive, printed element stress tables in the new version 
of SCARE are optional as shown in Fiq. 3. In addition, 
two new element cross-reference tables are printed. 

The first table lists the shell element number and 
gives the corresponding solid element with which its 
grid points are shared. The second table lists the 
solid element identification number and gives up to six 
associated shell elements (a HEXA element could have 
all of its six faces as an external surface). Element 
areas and temperatures (nodal averaged) are summarized 
in the following table. Next the selected fracture 
model is identified and the room temperature statisti- 
cal fracture parameters are shown. Additionally, a 
table of discrete input temperatures with corresponding 
material parameters, which were either internally cal- 
culated or directly supplied, is printed. For the 
shear-sensitive fracture models, the crack shape is 
identified along with a more specific description of 
the fracture criterion. The last table in the SCARE 
output file contains an element results summary, list- 
ing the shell element number, corresponding element 
survival and failure probabi 1 ities, and interpolated 
(based on temperature) element material parameters. 
Finally, the overall component probability of failure 
as well as the component probability of survival are 
printed. 

THEORY 

The most widely used mathematical models describ- 
ing the statistical nature of fracture in brittle mate- 
rials have been previously summarized by Gyekenyesi 
(1986). It is generally recognized that the weakest 
link hypothesis and the two-parameter Weibull cumula- 
tive strength distribution are appropriate to model inq 
surface flaw induced failure in ceramic structures. 
Consequently, the uniaxial fracture data is approxi- 
mated by 
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p 


fs ~ 



dA 


( 1 ) 


where the subscript s denotes parameters' associated 
with the surface, Pf s is the cumulative failure 
probability, a os is the scale parameter with 

1/m 

dimensions of stress x (area) , m s is the Weibull 


modulus which measures the degree of variability, a s 
is the applied surface tensile stress, and A the 
stressed area. The scale parameter is often called 
the unit area characteristic strength, and as m s 
increases it approaches the material's ultimate 
strength. The selectively used Weibull threshold 
stress parameter, o us , is taken to be zero in Eq. (1) 
and is not shown therein. 

In the analysis of failure of brittle materials 
subject to multiaxial stress states, the Weibull model, 
when combined with the PIA hypothesis, yields 



where a\ s and < 12 $ are the principal in-plane 
stresses acting on the surface of the structure. 
Equation (2) has been widely used in the past to esti- 
mate failure probabilities of ceramic components (Wertz 
and Heitman, 1980). The failure probability using the 
normal stress averaging method, as extended to surface 
distributed flaw problems, can be calculated from 



where k W p S is the polyaxial Weibull crack density 
coefficient for surface flaws. This constant can be 
obtained by making the result of integrating Eq. (3), 
using the normal stress a n$ distribution on an 
arbitrary plane, obtained from the reduced, plane- 
stress Cauchy infinitesimal tetrahedron in principal 
stress space as shown in Fig. 4, for uniaxial stress 
cases, agree with the results obtained from the uni- 
axial, two-parameter Weibull equation. Unlike for 
volume flaw based analysis where closed form integra- 
tion is possible, the polyaxial Weibull surface crack 
density coefficient must be evaluated numerically. The 
contour integration is performed on the circumference 
of the unit circle C where the normal stress is 
tensile and neglecting regions where the normal stress 
is compressive. The crack-like flaws can then be 
regarded as located in these arbitrary planes which are 
tangent to the unit circle and are acted upon by o ns 
which is induced by the principal stresses and 

o 2 S . Since Eq. (3) is just the shear-insensitive 
case of the more general Batdorf polyaxial stress frac- 
ture model, its SCARE implementation follows a somewhat 
different format. The polyaxial Weibull equation has 
also been selectively used in the past (Paluszny and 
Wu, 1977), but since it neglects the effects of shear 
loads, it also underestimates failure for the more 
general loading condition. These unconservative pre- 
dictions have been recognized early in ceramic struc- 
tures and the use of amplification factors was proposed 
to better correlate predictions with measurements. 


In the previously described two multidimensional 
stress fracture models, no mechanistic fracture crite- 
rion was needed to predict impending failure. In prior 
work by Batdorf (1973), Batdorf and Heinisch (1978), 
Shetty et al., (1984), and Samos (1982), attention is 
focused on cracks and their failure under stress. 

Since there is yet no consensus regarding how to treat 
mixed-mode fracture, even in metallic structures, the 
SCARE program includes several fracture criteria and 
flaw shapes. These criteria and flaw configurations 
include those traditionally treated, although recent 
research in ceramic fracture emphasizes more advanced 
fracture criteria and better modeling of flaw shapes 
(Petrovic, 1985). 

Consider now a small uniformly stressed material 
element of area aA. The probability of failure 
under an applied state of stress can be written as 


p fs = Pis p 2s 


(4) 


where P\ s is the probability of existence in aA 
of a crack having a critical stress in the range of 
°cr to a cr + ^ a cr* an ^ ?2s denotes the probability 
that a crack of critical stress a cr will be ori- 
ented in a direction such that an effective stress <? es 
equals or exceeds o cr . o cr is defined as the 
remote, uniaxial, normal fracture stress of a given 
crack. Failure will occur when the effective stress (a 
function of chosen crack configuration and fracture 
criterion) exceeds o cr for a particular crack. P| s 
has the form 


dN 

P ls = 4A dT 


dcr 


cr 


(5) 


and 


P 


_ 

2s “ it 


( 6 ) 


where N s (o cr ) is the surface crack density function 
(the density of cracks having a critical stress <o cr ) 
and til is the angle in principal stress space contain- 
ing all the crack orientations for which ap S > a cr . 
Using the weakest link theory, the overall faiTure 
probability can be calculated from (Samos, 1982) 


r c n* 



The surface crack density function N s (o cr ) is a 
material property which is expected to be of Weibull 
form since it is directly the cause of the observed 
Weibull strength distribution. However, for values 
of m s » 1.0, and assuming that large number of 

cracks are sampled, that is o cr << o os A”^ m s, the 
Weibull distribution can be approximated by a power 
function, leading to 


m 

N s («cr) = k Bs°cr 

The surface flaw distribution parameters kg s and 
m s can be evaluated from experimental data using 
uniaxial specimens and an appropriate fracture crite- 
rion. Since cracks in a real material are assumed to 
have an equal likelihood of orientation, and it is 
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believed that there are so many cracks present that 
there is always one crack with the least-favorable ori- 
entation and that crack growth from this flaw causes 
failure, mode I fracture is intrinsic to uniaxial load- 
ing and the critical normal stress fracture criterion 

^ 1 . . . a . J J L. ^ m ^ k at A f nr ^ r 

U I wajfo LU UCtCI 1 1 ^JVJ. 

However, finishing operations such as grinding and 
environmental damage from erosion or oxidation may lead 
to anisotropic surface conditions, which then must be 
carefully accounted for in material characterization. 

It is important to note that the use of various 
mixed-mode fracture criteria to evaluate from 

uniaxially loaded specimens will result in different 
corresponding values of this flaw density parameter, 
from which it follows that the component reliability 
predictions will also change (Shetty et al., 1984). In 
fact the use of an identical mixed-mode fracture crite- 
rion to evaluate kg s from MOR bars first, and then 
subsequently calculate multiaxially loaded component 
reliability with it will result in contradictory 
trends, where with increased shear-sensitivity, compo- 
nent failure rate predictions will decrease. In the 
present version of SCARE, failure rate predictions 
increase with increasing shear-sensitivity, which 
follows from always assuming mode I fracture in uni- 
axial tests. Limited experimental data with monolithic 
ceramics (Samos, 1982) tend to validate the present 
approach, and the need for more pessimistic failure 
predictions than those obtained from the early phenom- 
enological models is evident from component test 
records. 

The statistical analysis of fracture is greatly 
simplified by assuming that cracks are shear- 
insensitive. For this case fracture occurs when 
°es = °ns >. a cr and there is no need to specify the 
crack shape or the materials Poisson's ratio. The 
effective stress just equals the normal stress acting 
on the crack plane. 

When cracks are shear-sensitive, equations for 
o es using various fracture criteria can be derived by 
equating selected parameters, induced by a uniaxial 
normal stress acting on a specified crack shape, to the 
selected parameters calculated for the mixed-mode or 
multiaxially loaded identical crack configuration. 

Some of the most widely used parameters to predict 
impendinq failure include the total coplanar strain 
energy release rate or crack extension force G c , 
the maximum strain energy release rate G max , the 
minimum strain energy density S c , the normal stress 
o ns , and the maximum tangential or hoop stress a ee . 

The G max , S c , and o 00 fracture parameters predict 
out-of-plane crack extension, and consequently, are 
believed to be the most accurate for mixed-mode crack 
propagation. In any event, for polyaxial stress 
states, the effective stress o ?s is a function of 
both CT ns and t s where t s is the shear stress 
in the crack plane. Both o ns and t s , acting on a 
plane whose normal and the maximum principal stress 
direction form an angle a, can be calculated from con- 
sidering the equilibrium of forces in Fig. 4. Summing 
forces in the normal and tangentical directions gives, 
respecti vely 


2 • 2 

a = Oi COS a + a 0 „sm a 


ns ' u ls^ 


2 s 


= (°ls - °2s) cos ° sin 


(9) 


Using the G c criterion and the two selected crack 
configurations, the effective stress equations obtained 
are (Samos, 1982) 


a es = Y a ns + T s ( Griffith crack) 

( 10 ) 


n L ^ + 1 ( ^ — . — 1 ( Cri f f i t h nntrh 1 

es y ns i - v y i.izib j ' 

where v is Poisson's ratio and S s is the out-of- 
plane shear stress since t s does not cause crack 
growth in a Griffith notch. In general, the total 
coplanar strain energy release rate can be expressed in 
terms of stress intensity factors, summing damage in 
all appropriate crack growth modes. From equality of 
cross-shears, S s = t s and Eq. (9) can be used 
to evaluate o es for both cracks. 

The angle w depends on the fracture criterion 
selected, the assumed crack configuration and on the 
applied stress state. Closed form expressions for <o 
can be derived for some fracture criteria in uniaxial 
and equibiaxial stress states (Samos, 1982). Assuming 
a uniaxial stress o i s = a s , and the normal stress 
(shear-insensitive) fracture criterion, we obtain from 
Eg. (9) at ^ns ~ ^cr* ® — ®cr* ^nd as 2ot^^» that 



Note that when a shear-sensitive fracture criterion is 
used, the crack shape must also be specified. In 
general, for two-dimensional stress states, tu must be 
determined numerically. Using the shear-insensitive 
case as an example, we obtain at fracture 

9 

°es = a ns = °cr = a ls cos a cr 

+ a 2s sin 2 a cr (12) 

If we define K - ^2s/ CT 1 s and we note that 
a cr = u>/2, Eq. (12) "yields 



In order to avoid singular conditions (o| s = a^) 
in Eq. (13), the following constraints must be imposed 
in calculating P 2 S : 

if ( ° C r / °ls " K)/(l-K) < 0.0 then w = it 

if (o cr /oi s - K) / ( 1-K) > 1.0 then w = o 

(“ = ’ if °ls ^ °cr 

if K = 1.0 then < 

(a. = O if c ls < a cr 

A similar procedure can be followed for the shear- 
sensitive Griffith crack using both of Eq. (9) in 
Eq. (10). However, for the Griffith notch it is compu- 
tationally convenient to define $ = cos 2 a cr , so that 
the fracture condition a es ^ a cr can be written 
as (Samos, 1982) 


c l* 2 + c 2* + C3 = 0 


(14) 
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Sol ving for <t> gi ves 


where 


<t> = COS a 


_c 2 *' 


V C 2 - 4C 1 C 


M „v2 /, 0.795\ 

c i ■ o - K > y -ttt-j) 

C 2 = (1 - K) [ 2 K ♦ (1 - K)^iJ 


c 3 = K - 


\ / 

From Eq. (15), a cr = cos _1 J*. Consequently 

r ’ 1 


U> 2 2 

— = — OL = — COS 

TT TT Cr Tf 


-c z ± _ 4 Cl c 3 


2c, 


1 12 


(15) 


(16) 


(17) 


equation for p£c in Eq. (7), the failure probabil- 
ity usinq the Batdorf approach for surface flaws can be 
calculated. Unlike for volume flaw analysis where a 
double numerical integration is generally required, 
surface flaw analysis uses only a single integration 
for the critical stress. Assuming that the stress 
state is constant in any given shell element (SCARE1 
version) or subelement (SCARE2 version), the area inte- 
gration over the material surface becomes trivial. For 
component reliability analysis, the elemental survival 
probabilities are summed according to established prob- 
ability axioms, that is, we evaluate the product of all 
the individual, constant stress state elemental reli- 
abilities. For two mutually exclusive events, such as 
the calculation of Pf 5 and P ss , it is well known 
that for a given area increment their sum is always 
unity, where P ss denotes surface survival proba- 
bility. Finally, if a component includes 
two concurrent modes of possible failure, such as due 
to simultaneous internal or external flaws, its total 
reliability can be determined from the product of its 
individual cause reliabilities. 


Following an analogous procedure to the above in using 
fracture criteria and flaw shapes, similar expressions 
for terms of principal stresses and a cr can 

be developed. Results of this analysis are summarized 
in Table I and details are given by Samos (1982). In 
addition, special stress states such as the equibiaxial 
loading case where K - 1.0 and when Eq. (17) becomes 
singular, are treated. 

The equations in Table I apply only when o ns is 
tensile. When a compressive normal stress state 
exists, the angle where the normal stress changes from 
tensile to compressive is calculated. This transition 
angle divides all possible crack orientations into a 
tensile normal stress region and a compressive normal 
stress region, where each region has a different effec- 
tive stress equation. In the compressive stress region 
a ns is set equal to zero, and with friction between 
crack faces ignored, a es is a function of only the 
applied shear load. In each region, the angle at which 
°es = a gr 1S calculated. This angle further divides 
the region into areas where ag S > a cr and o^ s < o cr . 
If there is no such angle within the qiven region, then 
either there is no contribution to 10 from that reqion 
or the entire region is included. A calculation com- 
paring a es to a cr for a selected angle a 
identifies which subregion, or whether the whole 
region, contributes to w. After determining each w 
from the tensile and compressive regions separately, 
the total angle w is simply the sum of the two 
contributions. 

An alternate approach to calculating p£c has 
also been previously developed, where the angle a is 
incremented over the unit circle, and at each discrete 
point the contribution of o ns or a es is summed 
depending on whether they are equal or greater than 
a cr . Since o cr is always positive, compressive 
a ns or a es never contribute to fracture. 

The same method can also be used in volume flaw based 
failure analysis, where two angles must be incremented 
so that the stresses are checked at all points on the 
unit sphere. This procedure, when compared to the 
previous method used to calculate P 2 S , is computa- 
tionally much more intensive. Consequently, SCARE 
employs the more efficient approach and the equations 
summarized in Table I are coded to calculate critical 
crack orientations. 

In the SCARE program the required integration is 
performed by usinq Gaussian Quadratures. Employing the 
power function form of N s (n cr ) and the appropriate 


MATERIAL STRENGTH CHARACTERIZATION 

Uniaxial ceramic strength data is usually obtained 
from flexure testing simple geometry specimens, and 
recording the extreme fiber tensile stress in the most 
highly loaded area at fracture. The Weibull parameters 
a os and m s as well as the Batdorf flaw density 
paramenter kg 5 are evaluated in SCARE using four- 
point MOR data with known qe^metry. Fracture is some- 
times assumed to always occur between the symmetrically 
placed inner loads on the tensile side (Samos, 1982), 
but in SCARE an effective area is used which permits 
failure in any tensile loaded areal element of the bar 
(Govila, 1983). It should be noted that for uniaxial 
tensile loadinq the effective area A e reduces to 
A, the specimen gauge surface area. In all other types 
of loading, A e is a function of specimen geometry 
and material Weibull modulus. If we express the uni- 
axial flexure failure probability in terms of the maxi- 
mum extreme fiber fracture stress of or MOR usinq 
Weibull form, that is 



then after evaluating C s and m s by the least- 
squares method, a os and kg s can be calculated. 
Using the effective area for a rectangular beam in 
four-point bending (Govila, 1983) 



where is the length between the outer loads, I? 

is the length between the inner loads, w is the beam 
width, and h is the beam height, we can rewrite 
Eq. (1) as 

P fs . 1 - exp |- 

Comparing Eqs. (18) and (20) where in both equations 
the maximum extreme fiber stress is used as the refer- 
ence value, we conclude that for the same failure 
probabi 1 ity, 



6 



ORIGINAL PAGE IS 

OF POOR QUALITY 


o 


OS 



( 21 ) 


In addition to determining a os and m s , the 
‘sfARF nrnnram rpnnirps knowledge of ko^. Assumina 
shear-insensitive fracture, Eq. (11) can be 
substituted into Eq. (7) to obtain 


where A e b is the effective equibi axial l.y loaded 
tensile area and the equibiaxial failure rate is 
assumed to also follow Eq. (18). When usinq the least 
squares method with Eq. (18) to evaluate C s and 
m s , Pf s can be determined from ranked fracture stress 
data using the simple equation 


Kis 1 “ TT+T 


fs '°js 


(28) 



Taking the natural log of Eq. (22) and integrating by 
parts yields 



where A e is given by Eq. (19) (A e « wL^) for the 
four-point bend specimen. Equation (23) is an Abel 
integral equation with the solution (Batdorf and 
Chang, 1977) 


A e N s 




(24) 


where N is the total number of uniaxial or equibi- 
axial tests and j is the rank of the surface fracture 
stress a j s . 

EXAMPLE 

In order to test the surface crack distributed 
reliability analysis in SCARE, several example problems 
were analyzed from the open literature. In addition, 
the rotationally loaded, silicon nitride annular disk 
(Gyekenyesi, 1986) was reanalyzed assuminq that frac- 
tures in the MOR bars as well as the disks were caused 
by surface flaws. The Weibull modulus m s was set 
equal to m = 7.65 as before, but the previously used, 
volume flaw based scale parameter a 0 , was readjusted 
so that Pf - Pf s in the MOR bars, when both are 
loaded in uniform tension. As a result of this 
requirement, we obtain for the four-point bend speci- 
men that. 


w 


2(m s w + w + h) 


ll/m. 


(29) 


x V° cr -"“Is do ls 

From Eq. (18) with o^ s = a f as the reference 
stress 




m -2 

[ 1n P ss M ■ m s fm s - 

- C oi 
sis 


(25) 


Introducing a new variable n defined as 
n = (2oj s /o cr ) - 1 and substituting Eq. (25) into 
Eq. (24), we obtain 


3/2-m 


N s (o cr ) - 


C s m s (m s ‘ 


1)a c S r f 

*1 


x (1 


m -2 


i) s (1 - n)^^ dn 


(26) 


Now from Eq. (21) C s = A e a Q ^ , and consequently 

and A 0 can be eliminated from Eq. (26) and 

used instead. Since N s (a Cf ) = k gs a c ®, k B$ can be 

evaluated by numerically integrating Eq, (26) using 
Gaussian quadratures. Since a« s is used, the 
effective area given by Eq. (19) is inherent in N s (a cr ). 
It is interesting to note that Mo cr ) can be deter- 
mined more easily from equibiaxial tensile tests since 
for = ct£ s = o ns = a, uli t = 1.0 and 


(o cr)=-AT ln(> ss=A 


eb 


eb 


a cr k Bs °cr 


(27) 


From Gyekenyesi (1986), a 0 = 74.82 MPa m^-3922 
(45,800 psi in. 0.3922) an( j the "A-size" test bars had 
dimensions of w = 0.635 cm (0.25 in.), h = 0.3175 cm 
(0.125 in.), L\ = 1.905 cm (0.75 in.) and 
12 = 0.9525 cm (0.375 in.). A e is given by Eq. (19) 
and V e is the effective volume in four-point bending, 
that is 


V e = 


L 1 + m s L 2 


(m $ + ir 


(30) 


Using the given data in Eq. (29), we obtain that 
a QS = 232.0 MPam 0-2614 (87,891 psi in. °* 2614 ). The 
dimensions of the disk are given in Fig. 5, and the 
same MSC/NASTRAN analysis was performed as for volume 
flaws with the exception that shell elements were 
added to identify external surfaces, to calculate 
their areas and to identify appropriate surface stress 
states. Reliability calculations were made at various 
speeds using several surface crack fracture models. 
Selected results from these analyses are shown in 
Fiq. 5 and Table II. For a given speed, failure prob- 
abilities are considerably less than those obtained by 
Gyekenyesi (1986) for all fracture models, indicating 
that failure was most likely due to volume flaws as 
originally concluded. Obviously, the main reason for 
the greatly decreased failure estimates is the much 
higher equivalent surface Weibull scale parameter o 0S . 
Results for the shear-sensitive Griffith crack and 
notch were practically identical, with the notch beinq 
slightly more shear-sensiti ve as expected. The impor- 
tance of post-mortem fractography to identify the 
nature of the fracture causing flaws is evident from 
the two widely different set of answers obtained for 
the same problem. Furthermore, the difference between 
predictions from various available surface crack frac- 
ture theories for a specified loading is considerably 
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less than the difference found in volume flaw based 
fai lure models. 

In addition to comparing surface and volume flaw 
based analyses for the same ceramic structure, failure 
probability was calculated for a transversely loaded 
circular alumina disk and compared to data obtained by 
Brockenbrough et al. (1985). This data consisted of 
calculated failure probabilities obtained from using 
Batdorf‘s shear-i nseni ti ve fracture theory in conjunc- 
tion with the finite element analysis code ANSYS. The 
selected disk had an outside radius r 0 = 25 mm 
(0.984 in.), a thickness of t = 2.5 mm (0.098 in.), 
was loaded by a circular line load at r = 1.5 mm 
(0.059 in.) and was simply supported at a radius of 
22 mm (0.866 in.). After matching the stress solu- 
tions for a given load from the two different finite 
element codes (MSC/NASTRAN and ANSYS), the failure 
probabilities were calculated at various loads. 
Typically, to obtain a Pf s = 0.50, SCARE2 predicted 
a total line load of 1105 N (248.3 lb), while accord- 
ing to Brockenbrough et al . (1985) (corrected version) 
a line load of approximately 1130 N (253.9 lb) had to 
be applied, showing good agreement between the two load 
predictions. 

CONCLUSIONS 

The general purpose, statistical, fast fracture 
failure probability code SCARE has been enhanced to 
include the possibility of failure due to two con- 
current but not interacting flaw populations. In addi- 
tion, reliability analysis with linear finite elements 
is now also possible. The program includes a number 
of widely used polyaxial fracture models, appropriate 
extreme value statistics and the ability to calculate 
material failure distribution parameters for both vol- 
ume and surface distributed flaws. Current work 
includes the planned addition of more advanced increas- 
ingly shear-sensitive failure criteria which permit 
out-of-plane crack extension for both types of imper- 
fections. In addition, the influence of bimodal flaw 
populations on individual material parameters, through 
censored data analysis, needs to be determined. 

Finally, the problem of a transversely loaded circular 
plate will be anewly studied to resolve some of the 
contradictory trends reported by Rufin et al. (1984) 
and Shetty et al . (1984). 
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TABLE II. - EXAMPLE 1 FAILURE PROBABILITIES AS A FUNCTION OF ROTATIONAL 
SPEED FOR VARIOUS FRACTURE MODELS 


[m s = 7.65; o 0S = 232.0 MPa (m) 0 - 2614 ; k Bs oos = 4.986; NGP = 15] 


Angular 

speed, 

rpm 

SCARE 2 (HEXA AND QUAD8 elements) 

Experimental 
(Ford Motor Co. ) a 
(Gyekenyesi, 1986) 

Weibul 1 
PIA 

Batdorf shear- 
insensiti ve 

Batdorf shear- 
sensiti ve G.C. 
G c criterion 

Batdorf shear- 
sensiti ve G.N. 
G c criterion 

80 000 

0.0031 

0.0034 

0.0044 

0.0045 

0.2017 

85 000 

.0078 

.0087 

.0113 

.0114 

.3367 

93 000 

.0304 

.0339 

.0438 

.0441 

.6321 

100 000 

.0898 

.0996 

.1276 

.1283 

.8714 

105 000 

.1796 

.1979 

.2496 

.2508 

.9640 

110 000 

.3309 

.3610 

.4418 

.4436 

.9949 

115 000 

.5497 

.5889 

.6857 

.6877 

.9997 

120 000 

.7828 

.8176 

.8908 

.8922 

1.0000 

125 000 

.9417 

.9579 

.9838 

.9842 

1.0000 

130 000 

.9945 

.9970 

1.0000 

1.0000 

1.0000 


G.C. - Griffith crack 
G.N. - Griffith notch 

Experimental data fitted to a Weibull equation and calculated at the speeds shown. 



Figure 1. - Computational elements of the SCARE2 reliability analysis program. 









Figure 2. - Flowchart for subroutine ANGLE based on surface flaw fracture analysis. 
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IPRINT 

CONTROL INDEX FOR PRINTING STRESSES 
1: DO NOT PRINT ELEMENT STRESSES 
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END OF 

FIRST CARD OR STATEMENT NUMBER ONE 


BEGINNING OF SECOND CARD OR STATEMENT NUMBER TWO 

1-5 

Kan 

NUMBER OF TOTAL HEXA AND PENTA ELEMENTS 

6-10 

■m 

NUMBER OF HEXA ELEMENTS USED 

11-15 

NP 

NUMBER IF PENTA ELEMENTS USED 
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21-25 
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NT 
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TEMPERATURE 

41-45 
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ENTIRE STRUCTURE 
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NUMBER OF TEST TEMPERATURES AT WHICH 
MATERIAL DATA IS SPECIFIED 

56-60 

LONL 

CONTROL INDEX FOR SELECTING USE OF LINEAR 
OR NONLINEAR (QUADRATIC) ELEMENTS 
0: ONLY LINEAR ELEMENTS ARE USED IN THE 
MODEL 

1: ONLY QUADRATIC ELEMENTS ARE USED IN 
THE MODEL 


Figure 3. - SCARE surface flaw based analysis master control 

DECK DATA REQUIREMENTS. 














cos a 


m = sin a 

2 • 2 

o ns = oi s cos c a + o 2 S sin 4 - a 
T s = (Ois - 02 S )cos a sin a 


Figure 4. - Stresses on the reduced plane-stress Cauchy 

INFINITESIMAL TETRAHEDRON IN PRINCIPAL STRESS SPACE. 


PROBABILITY OF FAILURE 



Figure 5. - Example 1 probability of failure versus disk 

ROTATIONAL SPEED FOR VARIOUS FRACTURE MODELS (SCARE 2 

data). Surface flaw induced failure is assumed in all 
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